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The collimation efficiency for 208p^82+ 

ion beams in the LHC is predicted to be lower than 
requirements. Nuclear fragmentation and electromagnetic dissociation in the primary collimators 
create fragments with a wide range of Z/A ratios, which are not intercepted by the secondary 
collimators but lost where the dispersion has grown sufficiently large. In this article we present 
measurements and simulations of loss patterns generated by a prototype LHC collimator in the 
CERN SPS. Measurements were performed at two different energies and angles of the collimator. 
We also compare with proton loss maps and find a qualitative difference between 208pjj82+ j^^^^ 
and protons, with the maximum loss rate observed at different places in the ring. This behavior 
was predicted by simulations and provides a valuable benchmark of our understanding of ion beam 
losses caused by collimation. 

PACS numbers: 41.85. Si, 29.20.dk, 41.75.Ak 



I. INTRODUCTION 

The Large Hadron Collider (LHC) [Ij, presently being 
commissioned at CERN, will collide beams of protons 
and later heavy nuclei [7, starting with 208p|-j82+^ 
ergies never reached before. The main parameters of the 
beams of the two species are listed in Tab. |l] The LHC 
uses superconducting magnets, which operate with a high 
magnetic field at 8.33 T, near the quench limit, meaning 
that even a small temperature rise, of the order of a 1 K, 
can make the magnets pass from a superconducting to a 
resistive state. At the same time, with a stored proton 
beam energy of 362 MJ (see Tab.|l]), even a small beam 
loss of 4x 10'' protons in a magnetic element might induce 
enough heating to cause a quench. Therefore, all beam 
losses need to be tightly controlled and, for this purpose, 
a coUimation system has been designed [J [31 SJ 15] . This 
system is located in two warm insertions of the LHC and 
cleans halo particles from the beam in a controlled man- 
ner before they are lost elsewhere. 

The LHC collimation system is primarily optimized 
for proton operation but will be used also during ion 
runs. The collimation is however predicted to be less 
efficient for ions than for protons [5] , although the stored 
beam energy is almost a factor 100 lower, as can be seen 
in Tab. |l] In order to understand why this is the case, 
we give a very brief summary of the physics processes 
occurring when particles traverse the collimator material 
in Sec. |TT] and then of the functioning of the collimation 
system in Sec. |III| 

To better understand the LHC ion collimation, we have 
made an experiment on ion collimation in the CERN Su- 
per Proton Synchrotron (SPS), which is the main subject 
of this article. The experimental setup is presented in 
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Sec. llVl and the simulation methods in Sec. [V] In the re- 
maining sections, the measurement results are compared 
to simulations and we also make a brief comparison with 
the case of protons. 



II. PHYSICS OF HEAVY IONS IN 
COLLIMATORS 

In order to understand the LHC ion collimation ineffi- 
ciency, we give a brief overview of the physical processes 
occurring when heavy nuclei traverse the collimator ma- 
terial. 

A short review of the passage of charged particles 
through matter can be found in Ref. [7^ and an exten- 
sive theoretical treatment of ions in particular in Ref. [S]. 
Here we highlight two important processes: the en- 
ergy loss through ionization, which is described by the 
well-known Bcthc-Bloch formula, and the change of di- 
rection through many small-angle scattering events, so- 
called multiple Coulomb scattering (MCS). Angular de- 
viations can also be caused by nuclear elastic scattering, 
which is a significant effect for protons but negligible for 

These processes are present for all charged particles, 
while some others are peculiar to heavy ions. An imping- 
ing nucleus may lose one or several nucleons, in particular 
neutrons, through electromagnetic dissociation (EMD), 
which is a process with a logarithmic energy dependence 
taking place in ultraperipheral collisions. For a review of 
this field, see Refs. [311101 EI]- The nuclei may also split 
up in smaller fragments through nuclear inelastic reac- 
tions (see Ref. [12' for measurements of nuclear fragmen- 
tation of ^"^^Pb^^^ions at 158 GeV/nucleon in a carbon 
target, a situation close to the SPS collimation exper- 
iment described in this article). The cross section for 
nuclear inelastic interactions is only weakly dependent 
on energy and the created fragments have a wide range 
of masses. In Fig. [TJwe show the total cross sections for 
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FIG. 1; The cross sections for creating heavy isotopes (mass 
number A > 140 and charge Z > 60) by a 106.4 GeV/nucleon 
208p^S2+|^gg^j^ impacting on a carbon target. The cross sec- 
tions were simulated by FLUKA. 



the creation of different isotopes, the sum of the elec- 
tromagnetic and hadronic part, by 106.4 GeV/nucleon 
208p|-j82-i-jQj^g in a carbon target, as calculated by the 
Monte Carlo program FLUKA [HI [H US]. This cor- 
responds to the situation in the collimator in the SPS 
experiment. 

In each interaction, a number of secondary particles 
are created that constitute the hadronic shower [16] . The 
final energy deposition in the material is to a large extent 
due to these secondaries. 



III. LHC ION COLLIMATION 

The LHC coUimation system consists of primary and 
secondary collimators |,17 , both made of graphite jaws. 
Particles at large amplitudes hit first a primary collima- 
tor, which is shorter and the limiting aperture of the 
machine. There they interact in such a way that they 
are either absorbed or scattered back into the beam or 
to an even larger amplitude. The latter particles form 
a secondary halo, which is intercepted and absorbed by 
the secondary collimators, possibly several turns later. 
The condition on the angular kick Ax' given by the pri- 
mary collimator for a particle to be intercepted by the 
secondary collimator is [T7j 



Ax' > 



' {Ni-Nf)eN 



(1) 



Here ejv is the normalized emittance, 7 the usual Lorentz 
factor, (3x the usual optical function at the primary colli- 
mator and A^i , N2 are the transverse positions of the pri- 
mary and secondary collimators, given in the number of 
a (the standard deviation of an assumed Gaussian beam 
distribution) from the center of the vacuum chamber. 

We can now understand the inefficiency of the LHC 
ion coUimation by comparing the length of material that 
an ion would need to traverse in order to hit the sec- 
ondary collimator with the interaction length of nu- 
clear fragmentation processes. Taking standard LHC 



TABLE I: The parameters for the p+ and 208pjj82H 
beams at colhsion. 



LHC 



Particle 


P+ 


208p|^82+ 


Energy / nucleon 


7 TeV 


2.759 TeV 


Relativistic 7 


7461 


2963.5 


No. bunches 


2808 


592 


No. particles/bunch 


1.15x10'' 


7 X 10' 


Transverse normalized 






RMS emittances 


3.75 /im 


1.5 fim 


Stored beam energy 


362 MJ 


3.81 MJ 


Peak luminosity 


10^^ cm"^ S-' 


10^'^ cm-^ S-' 



2"8pb^2+parameters in Eq. ^ {Ni = 6, N2 = 7, e^v = 
1.5 fim, 7101 = 2963.5, Px ~ 135 m), we see that the 
angular kick that a ^^^Pb^^^ion would need to receive in 
the primary collimator in order to hit the secondary colli- 
mators is approximately 7 /irad. To have this RMS angle 
of the ion distribution exiting a collimator, the required 
length of the jaw would be 2 m if we use the Gaussian 
approximation of the Moliere formula ^ to calculate the 
RMS MCS angle Oq: 



13.6 MeV ^ 

vp V So 



1 + 0.038 In 



(2) 



Here Zq is the charge number of the incident ion, p its 
momentum, v its speed, s the distance traversed in the 
material and Sq the radiation length. 

However, the length of the primary collimators is only 
60 cm, and typical trajectories inside them much shorter, 
while the interaction lengths for nuclear inelastic inter- 
actions and EMD are around 2.5 cm and 19 cm respec- 
tively [6^. Since the transverse recoil of the heavy frag- 
ments is very small, it is clear that the heavy ions that are 
not absorbed in the primary collimator are likely to be 
fragmented without having obtained a sufficiently large 
angle to reach the secondary collimators. 

The created fragments (e.g. ^°''Pb, 203^j g^^^ others) 
have a different magnetic rigidity Bp (1+6) from the main 
beam (magnetic rigidity is defined as p/Ze = Bp for an 
ion with momentum p and charge Ze that would have a 
bending radius p in a magnetic field B). Therefore they 
are bent and focused differently. The fractional rigidity 
deviation 5 is given by 



6 = 



ZoA 
AoZ 



1 



Ap 

Po 



- 1, 



(3) 



where {Zq,Aq) are the charge and mass number of the 
original ion, {Z,A) of the fragment and Ap/po the frac- 
tional momentum deviation per nucleon of the fragment 
with respect to the main beam. These ions follow the lo- 
cally generated dispersion function from the collima- 
tor and may be lost downstream in the machine where 
the horizontal aperture A^ satisfies 



dx 6 > Ax 



(4) 




FIG. 2: A collimator jaw of the same type as installed in the 
SPS. The figure is taken from 22 . 

Here we have assumed that the dispersive contribution 
to the orbit is much larger than the betatronic part. Be- 
cause of the diflFerent 6 of the created ion species, they are 
lost at different locahzed spots, making the machine act 
as a spectrometer. Some impact locations could be out- 
side the warm coUimation insertion, meaning a potential 
risk of quenching superconducting magnets. 

To study the LHC ion coUimation inefficiency, a se- 
ries of simulation studies have been done [6j [18] . Since 
a large fraction of the systematic error in those simu- 
lations comes from the generation and tracking of the 
fragmented ions, an experiment on ion coUimation in the 
CERN Super Proton Synchrotron (SPS) has been per- 
formed and the results, presented in the following, are 
compared to simulations, not only in terms of loss lo- 
cations but with the goal of reproducing the absolute 
value of the losses measured by the beam loss monitors 
(BLMs) . We also make a brief comparison with the case 
of protons. 

IV. EXPERIMENTAL SETUP 

The SPS is a synchrotron of 6.9 km circumference that 
now serves as the last injector for the LHC [19] but also 
provides beam for various fixed-target experiments. It 
accelerates protons up to 450 GeV and ^"^Pb^^^nuclei 
up to 177 GeV/nucleon before extraction. 

A prototype of a secondary LHC collimator has been 
installed in the long straight section 5 (LSS5) [20l[2T]. It 
consists of a pair of 1 m long CFC graphite jaws, which 
can be moved independently to intercept the beam in the 
horizontal plane. A similar jaw is shown in Fig. J2] and 
the installation in the SPS tunnel is shown in Fig. [S]"^ The 
collimator is located in a position with small horizontal 
/3-function and dispersion in order not to introduce aper- 
ture bottlenecks. The optical functions and magnetic 
elements in the vicinity of the installation are shown in 
Fig. |4] the horizontal aperture in Fig. |5] and the vertical 
aperture in Fig. |6] 

Moving the collimator into the beam creates losses in 
the ring, since no second collimator is present to ab- 
sorb scattered particles. There is however no risk for 
quenches, since the SPS does not use superconducting 
magnets. The losses are recorded by 216 BLMs placed 
around the machine |23j . The BLMs are ionization cham- 
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FIG. 3: The installation of the prototype LHC collimator in 
the SPS tunnel. The figure is taken from [22] . 




FIG. 4: The /3-functions of the SPS downstream of the colli- 
mator and the locally generated dispersion dx from the colli- 
mator. We show also the magnetic elements (MB=main bend, 
QF=focussing quadrupole, QD=defocussing quadrupole). 
The BLMs and the collimator are indicated by vertical lines. 



bers mounted close to lattice quadrupoles (the inner part 
of a BLM and an example of the installation are shown 
in Fig.|7]) and filled with N2 gas at 1.1 bar. They consist 
of parallel aluminum plates, acting as anodes and cath- 
odes. The losses are read out and integrated over every 
machine cycle (18 s for ions). Fig.|4]shows the s- values of 
the four BLMs (called BL520, BL521, BL522 and BL523) 
immediately downstream of the collimator. BL520 and 
BL522 are mounted in the vertical plane below the beam- 
line, while BL521 and BL523 are in the horizontal plane 
on the inside of the ring. 

Beam loss data were collected during 
208pb^^+dedicated ion runs in late 2007 and the main 
beam parameters are shown in Tab. [IT] Data were taken 
both with a circulating beam at 106.4 GeV/nucleon 
and with a 5.9 GeV/nucleon beam intercepted by the 
collimator directly after injection. The corresponding 
momentum during the magnetic cycle of the SPS in 
both cases is shown in Fig. |8] In the case of the low 
energy beam injected on the collimator, the data were 
collected at the injection plateau during the first second 
of the cycle. Before the magnets started to ramp, the 
beam was dumped. 
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TABLE II: The parameters of the beams in the SPS used for the coUimation experiment. 



Particle 




208p|^82+ 


P 


Energy /nucleon (GcV) 


106.4 


5.9 


270 


Injected intensity (particles) 


6-7x10'' 


5-9x10' 


^Qi2_^Ql3 


Horizontal normalized emittance (1 cr) 


1 ^m 


1 fim 


2.7 fim 


Vertical normalized emittance (1 cr) 


1 fim 


1 fim 


4 fim 


Collimator steps 


0.1-1 mm 


fixed 


0.1-1 mm 


Collimator angles 


mrad, 2 mrad 


mrad 


mrad 




5250 527.5 5300 .5325 5350 5375 5400 




-0.01 
-0.02 
-0.03 
-0.04 
-0.05 
-0.10 
-0.15 
-0.20 
-0.30 
-0.40 
-0.50 




FIG. 5: Dispersive orbits of fragmented ions produced in one 
of the collimator jaws, shown together with the aperture. The 
collimator is located at s = 5222 m at the left edge. The 
vertical dashed lines indicate the locations of the four BLMs 
closest downstream. A large fraction of the total losses occur 
at the aperture limitation at s = 5277 m. 
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FIG. 6: Vertical betatron orbits coming out of the collimator, 
starting at y = 0, for different starting values of y' (vertical 
momentum normalized by longitudinal momentum). The col- 
limator is located at s = 5222 m at the left edge. The vertical 
dashed lines indicate the location of the four BLMs closest 
downstream. Particles with large vertical angles are lost close 
to BL520. 



FIG. 7; The installation of a BLM (BL521) in the SPS tunnel, 
shown together with a closeup of the inner assembly of the 
chamber. 



ate significant loss signals. We performed measurements 
using two different angles of the collimator with respect 
to the beam axis (0 mrad and 2 mrad) as schematically 
illustrated in Fig. [9] in order to vary the effective length 
travelled by particles in the collimator. 

The circulating beam current was measured by a Beam 
Current Transformer (BCT). These measurements were 
used later to normalize the BLM signals by the number 
of lost particles. The BCT has a resolution of a few 10* 
charges and the BLMs can detect around 10^ lost charges 
nearby. Since typical losses were a few 10^ charges, the 
combined uncertainty on the normalized signal is around 
10%. In the case where the beam was injected on the 
collimator, the normalization of the losses by the BCT 
could not be performed, since the injected intensity var- 
ied between different cycles. 

For comparison, data were also taken during high- 
intensity proton runs with the beam circulating at the 
flat top in the same way as for ^°*Pb*^^ions. These beam 
parameters are also presented in Tab. [IT] 



For the high energy measurements, the beam was in- 
stead accelerated up to 106.4 GeV/nucleon at the flat 
top, where the data-taking system continued to cycle 
while the magnets stayed at this energy and the RF sys- 
tem kept the beam at a constant energy. In this way the 
beam was circulating in the machine, while we scraped it 
with the collimator in several steps until all particles were 
lost. Because of the low intensity of the ^''*Pb*^~''beam, 
the collimator jaws had to be moved in close to the core of 
the beam (typical values are around 1 a) in order to cre- 



V. SIMULATION TOOLS 

In order to simulate the particle propagation through 
a lattice together with the particle-matter interactions 
in the collimators, a specialized program, ICOSIM (for 
Ion COUimation SIMulation), has been developed 
ICOSIM reads files from MAD-X [23 defining the op- 
tics and aperture of a machine, in order to make it 
straightforward to simulate any machine for which such 
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FIG. 8: The ^''^Pb^^+momentum per nucleon during three 
cycles of the SPS in the case of acceleration and a circulating 
beam at constant energy (CBCE). In the acceleration cycle, 
the dumping of the beam during the collimation measure- 
ments is indicated, as well as the extraction momentum. 




FIG. 9: Schematic sketch of collimator jaws kept parallel to 
the beam direction when moved in (left) and with a 2 mrad 
angle (right). 



a representation exists. The program creates an ini- 
tial beam distribution that is tracked in the variables 
{x,x' ,y,y' ,S, A, Z) through the lattice. The number of 
particles in the distribution is not constant, since inter- 
actions in the collimators may split particles into several 
smaller fragments with different values of A, Z and par- 
ticles hitting the physical aperture are removed from the 
tracking. Synchrotron oscillations are neglected, because 
typical time scales are long (Trf ~ 500 turns for the 
LHC and Tjip « 440 turns for the SPS) compared to the 
time between a first interaction of an ion with a collima- 
tor and the final loss (1-10 turns for most particles, 100 
turns in rare cases). 

Particles are tracked up to the first collimator impact 
using a linear matrix formalism between collimators with 
a small artificial emittance blowup on every turn, which 
accounts for diffusion under single beam effects, which 
need not be specified. The linear tracking is typically 
done for 10^ turns and serves to generate impact coordi- 
nates on the collimator. This method was developed with 
the aim of simulating fixed jaws and a growing envelope. 
We have chosen to use the same method to simulate the 
SPS experiments, although in this case the jaws were 
moving into the beam with a speed of 4 mm/s. This 
choice was made with the intention of benchmarking the 
specific generation of impact coordinates on the collima- 
tor and since the steps of the jaw were small compared 
to the discontinuities in the SPS aperture. As explained 
later, this fact makes the ratio between particles lost at 



different impact locations almost independent of the jaw 
positions. Thus, it is a very good approximation to use 
this simulation method also in the SPS. For the same rea- 
son, impact positions are relatively independent of orbit 
distortions and therefore we use the perfect machine op- 
tics in the tracking. The diffusion of the envelope for 
the SPS simulations was chosen to 92 nm/turn, corre- 
sponding to the jaw speed. The resulting average impact 
parameter (defined as the distance between particle im- 
pact and the inner edge of the collimator jaw) is 22 /xm. 

From the first collimator impact, element by element 
tracking is used in ICOSIM and chromatic effects at lead- 
ing order, as well as sextupoles in thin kick approxima- 
tion, are included. Higher order multipoles are neglected, 
since particles are typically lost in less than 100 turns, 
which is not enough to see a significant effect of small 
non-linearities. At the end of each element, a check is 
performed to determine which particles are outside the 
aperture. The impact momenta and coordinates on the 
vacuum chamber are estimated by a linear interpolation 
inside the magnetic element. 

Two different methods can be used to simulate the 
particle-matter interaction in the collimator: 

• FLUKA XS: ICOSIM has a simple built-in Monte- 
Carlo program, which includes multiple scatter- 
ing in a Gaussian approximation [7], ionization 
through the Bethe-Bloch formula, nuclear fragmen- 
tation and electromagnetic dissociation. The last 
two processes are simulated by sampling from tabu- 
lated cross sections, created by FLUKA for all pos- 
sible fragments and fragmentation channels. Only 
the heaviest fragment created in each interaction 
is tracked. In the LHC, only these fragments are 
important to follow, since ions with \6\ > 0.05 are 
already lost in the warm collimation insertion. 

• FLUKA full: The transport and interaction of all 
particles through the collimator region is evaluated 
on each turn by external calls to FLUKA, which 
simulates the full shower in a 3D model of the col- 
limator geometry. This method is more detailed 
and sophisticated but significantly slower in terms 
of computation time. 

The BLM signals depend not only on the number of 
ions lost nearby, but also on the mass of the ions (assum- 
ing the same energy per nucleon, a heavier ion carries 
more energy) , the distribution of impact parameters and 
the amount and type of material they have to traverse 
before reaching the monitor. At some BLMs, with less 
nearby material, particles lost far away may cause a sig- 
nal, while BLMs that are well shielded by magnetic ele- 
ments may only see small traces of the showers caused by 
the closest losses. In order to accurately simulate this for 
a quantitative comparison with data, the particle-matter 
interaction of the lost ions in the full geometry needs to 
be taken into account. 

As discussed later, the main loss location is right down- 
stream of the collimator. Thus, the 3D geometry of the 
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FIG. 10: The 3D geometry as implemented in FLUKA around 
a BLM in the vertical plane (BL520, top) and in the horizontal 

plane (BL521, bottom). A photo of the central part of the FIG. 11: The geometry as implemented in FLUKA around the 
real geometry around BL521 is shown in Fig. [7] monitor BL521 with the simulated energy deposition from a 

typical loss superimposed. The loss is indicated by the dashed 

arrow. 



TABLE III: The measured and simulated BLM signals (in 
Gy/10^° lost particles) during 106.4 GeV ^"^^Pb^^+operation, 
as shown in the top part of Fig. |13[ 



method/BLM 


520 


521 


522 


523 


meas. mrad 


0.079 


0.73 


0.01 


0.030 




±0.012 


±0.18 


±0.009 


±0.009 


FLUKA full mrad 


0.027 


0.52 


0.0008 


0.024 


FLUKA XS mrad 


0.0 


0.55 


0.0008 


0.032 


meas. -2 mrad 


0.079 


1.02 


0.007 


0.034 




±0.009 


±0.15 


±0.009 


±0.0019 


FLUKA full mrad 


0.031 


0.84 


0.0007 


0.026 


FLUKA XS mrad 


0.0 


0.73 


0.001 


0.036 



magnetic elements around the monitors BL520, BL521, 
BL522 and BL523 was implemented in FLUKA, as illus- 
trated for BL520 and BL521 in Fig. [lO| The magnetic 
field in the magnets nearby was neglected. A simulation 
including the magnetic field of the quadrupole nearest to 
BL520 showed a negligible difference with respect to the 
case without any field. The momenta and impact coor- 
dinates on the inside of the vacuum pipe of all particles 
lost within a 15 m interval of each BLM were recorded 
in ICOSIM and fed as starting conditions into FLUKA 
and the resulting energy deposition in the N2 gas inside 
the BLMs was converted to dose in Gy. An example of 
the simulation of the shower caused by the lost particles 
close to BL521 is shown in Fig. [TT] 

In the ICOSIM simulations of the SPS, typically more 
than 2 x 10^ particles were tracked in order to arrive at 
a statistical uncertainty of 1-5% on the final simulated 
BLM signal, depending on BLM. 



VI. RESULTS: 106.4 GEV/NUCLEON, 
PARALLEL JAWS 

We focus first on the case of a 106.4 GeV/nucleon cir- 
culating ^"^Pb^^'^beam and parallel collimator jaws and 
compare with other cases in later sections. 

A typical example of a loss map measured during a 
machine cycle with these conditions is shown in Fig. 12 
together with the corresponding simulated loss map from 
the FLUKA full method. The detector background, con- 
sisting of noise and other beam losses not caused by the 
collimator movement, had to be subtracted. As back- 
ground we used the loss map from the machine cycle 
before the collimator movement. A similar approach was 
already used in Ref. [21 to benchmark the SixTrack pro- 
gram p5l for proton beams. The only stable loss lo- 
cations, clearly separable from the background, are just 
downstream of the collimator, both in simulations and 
measurements. This holds true also with jaws tilted by 
2 mrad and at 5.9 GeV/nucleon. Therefore, in the re- 
mainder of this text we focus on the four BLMs down- 
stream of the collimator which see the highest signal. 

The top part of Fig. [13] and Tab. |III| show the average 
measured BLM signals, normalized to 10^° lost particles 
(using the BCT) and averaged over several machine cy- 
cles, together with the corresponding simulation results 
for the different methods of representing the particle- 
matter interactions in the collimator, which we discuss 
in Sec. |IX[ The ion loss pattern from the measurements 
is qualitatively very similar to the simulations, with the 
maximum signal on BL521. 

The fact that almost all losses occur close to the col- 
limator, as well as the ratio between signals on different 
monitors, can be qualitatively understood by consider- 
ing first the simulated (5-spectrum of nuclei leaving the 



collimator, which is shown in Fig. 14 We show the num 
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FIG. 12: Example of simulated ICOSIM (top) and mea- 
sured (middle) loss map of the whole SPS ring with a 
106.4 GeV/nucleon circulating ^°*Pb*^+beam. The bottom 
part shows a closeup of the loss peak in the measurements, 
with the names of the BLMs with the highest signals indi- 
cated. The simulated losses were binned in 5 m intervals. 
The collimator is located at s = 5222 m, just upstream of the 
large loss peak. 
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FIG. 13: Average loss map for 106.4 GeV/nucleon 
^"**Pb*^+ions (top) and 270 GeV protons (bottom) with par- 
allel jaws (blue bars) and jaws tilted by 2 mrad (red bars). 
The measurements are compared with simulations using both 
methods of simulating the particle-matter interaction in the 
collimator, as described in text. All results were normalized to 
10^" lost particles. The standard deviation between different 
measurements is indicated. 



ber of nucleons originating from ions with a certain S, 
since the FLUKA simulations show that an ion imping- 
ing in a magnetic element is fully fragmented within a 
few decimeters, meaning that the BLMs intercept mainly 
secondary particles originating from the hadronic shower, 
which to good approximation is similar to the shower cre- 
ated by the corresponding number of free nucleons . 
Therefore the BLM signal caused by an ion lost in a spe- 
cific location is approximately proportional to its number 
of nucleons. The heights of the peaks in the figures give 
thus an approximate estimate of the fraction of the BLM 
signal caused by different values of S. 

In order to explain the loss pattern, we consider this 



i5-spectrum in Fig. 14 together with the dispersive orbits 
of ion fragments with these (5-values starting at one of the 
collimator jaws, shown in Fig.[5j Here several trajectories 
for typical values of 6 are presented, together with the s- 
values of the BLMs. Fragments with \S\ < 0.013 stay 
inside the vacuum chamber and can make a full turn in 
the machine, while particles with larger \6\ are lost deter- 
ministically downstream of the collimator. Fig. [5] clearly 
demonstrates the spectrometer effect discussed earlier, 
since fragments with different values of 6 are lost at differ- 
ent longitudinal positions. Values of 5 for some common 



fragments are shown in Tab. IV Therefore, each of the 



four BLMs considered in detail sees a different spectrum 
of ions, shown in Fig. |15| This can help us to understand 
the origin of the signals at each BLM: 

• BL520: The main loss mechanism in the vicinity 
is not dispersion but instead the vertical aperture, 
which intercepts particles with large betatron an- 
gles. This can be understood from the spectrum of 
vertical angles of the fragments exiting the collima- 
tor, as shown in Fig. [16) and typical vertical orbits 
of particles with large scattering angles as shown 
in Fig. [6] The horizontal aperture is significantly 
larger (see Figs, [s] and |6|, making the vertical aper- 
ture the limitation for large-angle particles. Mainly 
light fragments are lost close to BL520, consistent 
with the expectation that they receive significantly 
larger transverse recoils in the fragmentation pro- 
cesses. The expected signal is significantly lower 
than measurements, maybe because BL520 is lo- 
cated only 30 m downstream of the collimator. It 
might see traces of shower particles from the colli- 
mator, in particular high energy neutrons created 
by electromagnetic dissociation, for which we have 
not attempted a detailed modelling. We have es- 
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TABLE IV: The 5 from fragmentation of the most common 
ions created in the coUimator as calculated with Eq. [S] for 
Ap/po = 0. 



207p^ 


207rpj 


20«Pb 


206rpj 


205 


205 rpj 


-0.0048 


0.0075 


-0.0096 


0.0026 


-0.014 


-0.0023 




*He 


^He 


«H 






0.0053 


-0.21 


-0.41 


0.18 


-0.21 


-0.61 



timated this contribution through a FLUKA simu- 
lation of the neutron propagation through a rough 
model of the tunnel and around 45% of the sig- 
nal shown in Fig. 13 (top) on BL520 comes from 
showers induced by scattering neutrons. However, 
since the neutron contribution is very sensitive to 
objects in the tunnel acting as scatterers, which are 
not properly taken into account in the simulation, 
this should be considered as a rough estimate. 



BL521: Ions with -0.2 < 5 < -0.08 are lost near 
the aperture limitation at s = 5277 m, close to 
BL521. Fig. [T4] demonstrates that this corresponds 



15 



to a large fraction of the fragments and Fig 
shows that BL521 sees a very wide spectrum of 
ions. Since BL521 is located only 2 m downstream 
of this position in negative x with almost no shield- 
ing material in between (see Figs, [t] [Toj and 111, 
this monitor is expected to show a very high sig- 
nal when ion beams are scraped with the collima- 
tor. Consequently, BL521 has the maximum signal 
both in measurement and simulations, which pre- 
dict a value around 30% lower than measurements. 



• BL522: BL522 is placed in the vertical plane, while 
the dispersive losses of mainly heavy fragments oc- 
cur in the horizontal plane. This makes the uncer- 
tainty of the shower simulation much higher, since 
the signal is caused by secondary particles very far 
away from the central core of the shower. BL522 
showed a much higher signal in measurements than 
simulations, although fluctuations between differ- 
ent cycles were very large (of the same order of 
magnitude as the measured signal), a fact that is 
not well understood. BL522 had the lowest signal 
(around 1.4% of the signal on BL521) of the BLMs 
that we consider in detail. 

• BL523: Only particles with a magnetic rigidity 
{5 « -0.02) close to the original 208PI382+ -^^ 
lost in the vicinity. Therefore the spectrum of lost 
ions consists mainly of heavy fragments. This sit- 
uation is similar to what can be expected in the 
cold regions of the LHC, since all particles having 
\6\ > 0.05 are already lost in the warm insertion. 
The signal from the FL UK A full simulation is 17% 
lower than measurements. 

It is clear from Fig. [5] that the expected distribution of 
losses is relatively independent of closed orbit distortions 
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FIG. 14: The simulated distribution of 5 calculated accord- 
ing to Eq. (|3|, using the FLUKA full method, of all high 
energy nuclei except ■^''**Pb*^"^ coming out of the collimator 
(top) and of the ions lost within a 15 m interval upstream of 
each BLM (bottom) at 106.4 GeV/nucleon with parallel jaws. 
The height of the bars show the number of nucleons belonging 
to ions having 5 within a certain interval. 



(which were not accounted for in the simulations) and the 
position of the collimator. Displacing the trajectories in 
the figure vertically by a few mm does not significantly 
change their impact positions in s, given the large impact 
angles and the aperture of the SPS. Likewise, the disper- 
sive orbits starting at the other jaw, located in negative 
X, have a similar longitudinal impact distribution. Thus, 
only the magnitudes of the measured signals, not their 
relative ratios, change when the collimator is moved in 
closer to the core of the beam and therefore intercepts 
more particles. This behavior was predicted by simula- 
tions and confirmed by measurements. 

The magnitudes of the simulated signals agree well 
with measurements (see Fig. 13 1 , although they are lower. 



To estimate the expected error, we assume approximately 
a factor 2-3 uncertainty of the FLUKA simulation of 
the energy deposition in the ionization chambers, since 
we consider secondary particles far from the core of a 
shower caused by ions with a grazing incidence. This is 
consistent with other comparisons between simulations 
and measurements of beam loss induced BLM signals at 
RHIC [27] and HERA [28]. 

Furthermore, the fragmentation cross sections in the 
collimator have significant uncertainties. A compari- 
son between simulation codes done in a separate study, 
which we intend to publish elsewhere, suggests that the 
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FIG. 15: The A distribution of ion fragments lost within 15 m 
upstream of each BLM at 106.4 GeV/nucleon and parallel 
collimator jaws as simulated with the FLUKA full method. 
The heights of the bars show the number of nucleons from 
ions having a certain A. 
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FIG. 16: The simulated distribution from the FLUKA full 
method of y' of all ion fragments coming out of the collimator 
(top) and of the ions lost within a 15 m interval upstream of 
each BLM (bottom) . The height of the bars show the number 
of nucleons belonging to ions within a certain interval. 



uncertainty of the largest cross sections to create spe- 
cific isotopes may well be 50%. Therefore, a new set 
of ICOSIM simulations were performed (using FLUKA 
XS), where we in each run resampled all cross sections 
with a 50% standard deviation around their initial value. 
This caused a standard deviation of 9% of the signal on 
BL521 and approximately 20% at the other BLMs, which 
see a narrower i5-spectrum and therefore are more sensi- 
tive to variations in single cross sections. 

There is also an uncertainty on the measured jaw 
angles, which has been estimated to be less than 
0.2 mrad [29 . Further FLUKA XS simulations show that 
this causes a 17% deviation on the BLM signal. Other 
systematic errors in the measurement were estimated at 
10% (see Sec.|lVl. 

However, these errors give only small corrections to the 
uncertainty of the FLUKA simulations. The same holds 
true for the single pass tracking from the collimator to the 
loss point, the influence of a non-perfect machine optics 
(as explained in Sec. V]) and for the impact coordinates on 
the collimator, which arc much better known in the SPS 
than in the LHC. The SPS collimator was moved into the 
beam close to the core, where the distribution is Gaussian 
with good approximation, while the LHC collimators will 
intercept halo particles, scattered to high amplitudes by 
beam dynamics processes that are hard to quantify. 

Altogether, we expect the simulation to be accurate 
within a factor 3 and conclude that the discrepancies be- 
tween measurements and simulations are within expected 
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FIG. 17: Measured and simulated ^°**Pb*2+loss maps with 
parallel jaws for a 5.9 GeV/nucleon beam injected on the 
collimator (green bars) and circulating 106.4 GeV/nucleon 
(blue bars). 



error margins, except for the highly fluctuating BL522. 
This error margin might be too pessimistic when con- 
sidering the BLMs in the horizontal plane (BL521 and 
BL523). 
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FIG. 18: The simulated A distribution of ion fragments from 
the FLUKA full method lost within 15 m upstream of each 
BLM at 5.9 GeV/nucleon. The heights of the bars show the 
number of nucleons from ions having a certain A. 



VII. 



RESULTS: 106.4 GEV/NUCLEON, TILTED 
JAWS 



Fig. 13 shows the average measured and simulated loss 
maps with the collimator jaws tilted by 2 mrad when 
moved into the beam (see Fig. [9]) . Compared to the case 
with parallel jaws the signal on BL521 increases by 40% 
in measurement, which is well reproduced by the simula- 
tions (34% in the FLUKA XS simulation and 60% in the 
FLUKA full simulation). This increase can be under- 
stood from the fact that the average distance travelled 
by the particles inside the collimator decreases, meaning 
that less particles are absorbed. 

Furthermore, the signal at BL523 shows a correspond- 
ing increase with tilted jaws both in measurements and 
simulations, while observed differences in the very small 
signal at BL522 are within the statistical error mar- 
gin. At BL520 an increase by 18% was predicted by 
the FLUKA full simulation but not reproduced by the 
measurements. 



VIII. RESULTS AT 5.9 GEV/NUCLEON 

In Fig. [17] we show the measured and simulated loss 
maps using a 5.9 GeV/nucleon ^''*Pb®^'''beam injected 
on the collimator with parallel jaws, together with the re- 
sults at 106.4 GeV/nucleon with parallel jaws presented 
earlier in Fig. |13[ Since the normalization to the inten- 
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FIG. 19: The simulated distribution of 5 calculated accord- 
ing to Eq. (|3| from the FLUKA full method of all ion 
fragments except ^"^Pb*'^"'" coming out of the collimator at 
5.9 GeV/nucleon. The height of the bars show the number of 
nucleons belonging to ions having S within a certain interval. 



sity decay measured by the BCT was not possible (see 
Sec. 



IV I , the normalization was instead done with re- 



spect to the sum of the four BLMs considered for both 
data sets and simulations. We noted however that the 
measured BLM signals dropped by around 5 orders of 
magnitude compared to the higher energy. 

As shown in Fig. [17] the fraction of the total signal at 
BL521 decreases at the lower energy, while it increases on 
BL520 and BL523. This is well reproduced by both sim- 
ulation methods, although they consequently predict too 
low a fraction on the BLMs in the vertical plane (BL520 
and BL522). The FLUKA full method also shows lower 
magnitudes of the signals than FL UK A XS. 

The change in loss pattern can be qualitatively under- 
stood from the energy loss in the collimator. Using the 
Bethe-Bloch formula, we see that a 106.4 GeV/nucleon 
208pj^82+|Q^ loses around 0.16 GeV/cm in carbon due 
to ionization, while a 5.9 A GeV ^"^Pb^^^^'^ion loses 
0.13 GeV/cm. This corresponds however to very differ- 
ent fractional energy losses: 0.15%/cm of the total energy 
for 106.4 A GeV ions and 1.9%/cm at 5.9 A GeV. Since 
typical paths inside a collimator are several centimeters, 
the energy loss caused by ionization becomes significant 
at low energy (see the dispersive orbits in Fig.[5|, mean- 
ing that the dispersion of the ions exiting the collimator 
is not only due to a different magnetic rigidity, but also 
because of a different energy per nucleon. The conse- 
quence is that most ions receive larger values oi \5\ than 
at the higher energy and are lost further upstream in the 
ring. Furthermore, ^°^Pb^^^ions that pass the collima- 
tor without fragmenting may receive a large enough \5\io 
make them lost close to all considered BLMs. This can 
be seen in Fig. [18] showing the mass spectrum of ions 
lost near each BLM. 

At the same time the production of fragments changes 
slightly: At 5.9 GeV/nucleon, the cross section for elec- 
tromagnetic dissociation goes down significantly |10j . 
while the cross section for nuclear inelastic interaction, 
on the other hand, is only weakly dependent on the en- 
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ergy [301 E]- Because of the different production rates 
and the ionization energy loss, the (5-spectrum of the ions 
leaving the collimator is different at 5.9 GeV/nucleon. 
This is shown in Fig. [T9| 



IX. COMPARISON OF SIMULATION 
METHODS 

In the FLUKA XS simulations no signal is seen at 
BL520 with the 106.4 GeV/nucleon beam, since only the 
heaviest fragment produced in each collimator interac- 
tion is tracked, while the signal at BL520 is caused by 
very light fragments (see Fig. 15). During operation at 



5.9 GeV/nucleon, heavy fragments are lost also at BL520 
(see Fig. 18 1, and consequently a non-zero signal is sim- 
ulated by the FLUKA XS method. 

Furthermore, the FLUKA full simulation at 
106.4 GeV/nucleon predicts some smaller loss peaks 
between BL520 and BL521, caused by very light parti- 
cles with a positive energy offset. These peaks are not 
present in the FLUKA XS simulation but are also too 
far from the BLMs to be seen in the measurements. 

However, at BL521, BL522 and BL523, where medium 
and heavy fragments are important, the FLUKA XS 
method gives an excellent agreement both for parallel 
and tilted jaws, while the gain in tracking speed is on 
the order of a factor 10 compared to the FLUKA full 
method. Furthermore, FLUKA XS reproduces well the 
qualitative changes of the loss pattern when the beam 
energy is changed. Therefore, we conclude that FLUKA 
XS is the preferred method for large machines with many 
collimators where light fragments are of low importance. 
This is the case for the LHC, since almost all light frag- 
ments are already lost in the warm regions. 



X. COMPARISON WITH PROTONS 

For the sake of comparison, proton runs in the SPS 
were also simulated with ICOSIM, using the FLUKA full 
method, and compared with measurements from Septem- 
ber 2007 with a coasting high-intensity beam (parameters 
are given in Tab. [ll|. More details and analyzes of these 
measurements will be published elsewhere. 

Average measured and simulated loss map are shown 
The ratio between the higher peaks agrees 



in Fig. 13 



well with measurements, although there is a factor 2 dis- 
crepancy in magnitude. This could, just as for ions, be 
explained with systematic errors dominated by the un- 
certainty in the shower simulation. Furthermore, the in- 
tensity was more than 4 orders of magnitude higher than 
in the ion runs, which might introduce additional losses 
not caused by the collimator. This is however not under- 
stood in detail. 

It is clear from Fig. [l3] that there is a significant qual- 
itative difference between proton and ion loss patterns: 
the maximum signal for protons was found on BL520 
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FIG. 20: The simulated distribution of S (top) and y' (bot- 
tom) of all protons coming out of the collimator. 



(closest to the collimator), while in the ion runs it was 
found on BL521. This can be understood from the fact 
that the S of the protons is much lower, since they cannot 
fragment, which can be seen from Fig. 



20 This means 



that large betatron angles caused by multiple scattering 
are the main loss mechanism instead of dispersion, al- 
though the spectrum of vertical angles is also more nar- 
row than for ions (see Fig. 20 1. The difference in loss 



pattern is a striking parallel to the expected behavior in 
the LHC, and the ability of the simulations to quantita- 
tively predict this behavior in the SPS provides a very 
valuable benchmark of the simulation studies which have 
been performed for the LHC. 



XI. CONCLUSIONS 

We have done measurements and simulations of 
^"^Pb^^^ion loss patterns induced by a collimator in the 
CERN SPS, in order to benchmark simulation tools used 
for the LHC. Qualitatively, the features of the loss pat- 
tern can be well understood from the particle-matter in- 
teractions in the collimator and the behavior of the dis- 
persion function downstream of it. A wide range of nuclei 
are created in the collimator due to fragmentation of the 
original ^°®Pb^^'''beam. The created isotopes have differ- 
ent Z/A ratios and follow different dispersive orbits until 
they are lost, making the machine act as a spectrometer. 

Quantitatively, the simulated loss distribution corre- 
sponds well to measurements. In terms of absolute BLM 
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signals, predicted and measured values agree to about 
30% for the largest losses, while discrepancies can reach 
a factor 3 for lower loss peaks (excluding the highly fluc- 
tuating BL522). This is consistent with expected un- 
certainties. We performed measurements with differ- 
ent beam energies and collimator angles and observed 
changes in the loss pattern, which were reproduced by 
our two-stage simulations with ICOSIM and FLUKA. We 
used two methods to simulate the particle-matter inter- 
action in the collimator — a fast simplified Monte Carlo 
and a full FLUKA shower simulation — and found a good 
agreement with measurements for both methods when 
the losses are dominated by heavy fragments close to the 
original ion. The simplified model is therefore well suited 
to predict limiting losses in the super-conducting mag- 
nets of the LHC with good accuracy. Losses consisting 
of light ions are on the other hand not treated by the 
simplified method, and if such losses are important the 
full shower simulation should be used. 

Furthermore, we have studied beam loss data taken 
with a proton beam and found that the loss patterns 
induced by the collimator in the SPS are qualitatively 



different for ^"^Pb ions and protons. This difference is 
well reproduced by simulations. The ions are lost mainly 
due to a change in magnetic rigidity caused by fragmenta- 
tion, while the protons are lost because of large betatron 
angles from scattering. 

These results confirm and strengthen our understand- 
ing of ion beam losses related to coUimation and are a 
vital test of our ability to make predictions for the LHC. 



XII. ACKNOWLEDGEMENTS 

We would like to thank A. Ferrari, M. Magistris, 
G.I. Smirnov and V. Vlachoudis for inspiring discussions 
and help with FLUKA. Furthermore, we are thankful to 
M. Stockner for discussions and for providing the techni- 
cal specifications and characteristics of the BLMs, and to 
G. Arduini, B. Dehning, M. Jonker, D. Kramer, D. Man- 
glunki, C. Zamantzas and the AB/OP group for their 
help during the measurements. We thank also L. Jensen 
for providing specifications of the instrumentation. 



[1] O.S. Briining, P. Collier, P. Lebrun, S. Myers, R. Ostojic, 
J. Poole, P. Proudlock et al. (Eds.), CERN-2004-003-V1 

(2004) . 

[2] J. Jowett, J. Phys. G: Nucl. Part. Phys. 35, 104028 
(2008). 

[3] R.W. Assmann et al, Proc. of the European Particle Ac- 
celerator Conf. 2004, Lucerne, Switzerland p. 536 (2004). 
[4] R.W. Assmann et al, LHC Project Report 850, CERN 

(2005) . 

[5] R.W. Assmann et al, Proc. of the European Particle Ac- 
celerator Conf. 2006, Edinburgh, Scotland p. 986 (2006). 

[6] H.H. Braun, Proc. of the European Particle Accelerator 
Conf. 2004, Lucerne, Switzerland p. 551 (2004). 

[7] S. Eidelman et al, Physics Letters B 592, 1 (2004), URL 
http : //pdg. Ibl . gov 

[8] J. Eichler, Physics Reports 193, 165 (1990). 

[9] G. Baur, K. Hencken and D. Trautmann, J. Phys. G: 
Nucl. Part. Phys. 24, 1657 (1998). 
[10] C.A. Bertulani, Int. Journal of Modern Physics 18, 685 
(2002). 

[11] C.A. Bertulani, S.R. Klein, J. Nystrand, Ann. Rev. Nucl. 

Part. Sci. 55, 271 (2005). 
[12] C. Scheidenberger et al, Phys. Rev. C 70, 014902 (2004). 
[13] A. Fasso, A. Ferrari, J. Ranft, and PR. Sala, CERN 

Report CERN-2005-10 (2005). 
[14] A. Fasso et al, Proc. of the Computing in High Energy 

and Nuclear Physics 2003 Conf., La JoUa (2003). 
[15] S. Roesler, R. Engel, J. Ranft, SLAG Report SLAC-PUB- 

8740 (2000). 

[16] A. Ferrari and P.R. Sala, Workshop on Nuclear Reac- 
tion Data and Nuclear Reactors, Miramare- Trieste, Italy 
(1996). 

[17] J. B. Jeanneret, Phys. Rev. ST Accel. Beams 1, 081001 
(1998). 

[18] R. Bruce, G. Bellodi, H.H Braun, S. Gilardoni and 



J.M. Jowett, Proc. of the European Particle Accelerator 
Conf. 2006, Edinburgh, Scotland p. 544 (2006). 

[19] M. Benedikt, P. Colher, V. Mertens, J. Poole, K. Schindl 
et al. (Eds.), CERN-2004-003-V3 (2004). 

[20] A. Bertarelh, O. Aberle, R.W. Assmann, S. Cala- 
troni, A. Dallocchio, T. Kurtyka, M. Mayer, R. Perret, 
S. Redaelli, G. Robert-Demolaize, Proc. of the Particle 
Accelerator Conf. 2005, Knoxville p. 913 (2005). 

[21] S. Redaelh, G. Arduini, R.W. Assmann, G. Robert- 
Demolaize,, Proc. of the European Particle Accelerator 
Conf. 2006, Edinburg h, Scotland p. 1810 (2006). 

[22] http : //Ihc-collimat ion-project .web. cern. ch/ 

[23] B. Moy, C. Ran, N. Aguilar, CERN LABII-RA/Note/73- 
9 (1973). 

[24] http://cern.ch/mad/', 

[25] F. Schmidt, CERN/SL/94-56-AP (1994). 

[26] S. Agosteo, C. Birattarib, A. Foglio Para, M. Silari and 
L. Ulrici, Nucl. Inst. & Methods A 459, 58 (2001). 

[27] R. Bruce, J. M. Jowett, S. Gilardoni, A. Drees, W. Fis- 
cher, S. Tepikian, and S. R. Klein, Phys. Rev. Letters 99, 
144801 (2007). 

[28] E. Holzer, B. Dehning, C. Fabjan, D. Kramer, M. Sapin- 
ski, and M. Stockner, Nuclear Science Symposium Con- 
ference Record, 2007. NSS '07. IEEE 2, 1679 (2007). 

[29] C. Bracco, Ph.D. thesis, Ecole Polytechnique Federale de 
Lausanne (2008). 

[30] K. Siimmerer, W. Briichle, D.J. Morrissey, M. Schadel, 
B. Szweryn and Yang Weifan, Phys. Rev. C 42, 2546 
(1990). 

[31] K. Siimmerer and B. Blank, Phys. Rev. C 61, 034607 
(2000). 



